2D bands and electron-phonon interactions in polyacene plastic transistors 
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We present a simple tight-binding model for the two-dimensional energy bands of polyacene field- 
effect transistors and for the coupling of these bands to lattice vibrations of their host molecular 
crystal. We argue that the strongest electron-phonon interactions in these systems originate from 
the dependence of inter-molecule hopping amplitudes on collective molecular motion, and introduce 
a generalized Su-Schrieffer-Heeger model that accounts for all vibrations and is parameter-free once 
the band mass has been specified. We compute the electron-phonon spectral function c? F{ijj) as 
a function of two-dimensional hole density, and are able to explain the onset of superconductivity 
near 2D carrier density n2D ~ lO^'^cm"^, discovered in recent experiments by Schon et al. 
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Recent studies |^J^ of high- mobility organic field-effect 
transistors by Batlogg and collaborators have demon- 
strated that their electronic properties can be tuned 
through an astonishing range, simply by adjusting the 
two-dimensional (2D) electron density n2D with a gate 
voltage. The high quality of these organic single crys- 
tals and of their interfaces with AIO3 dielectrics opens 
up new possibilities for studying the physics of coherent 
hand quasiparticlc transport in organic semiconductors. 
To date, however, the analysis of organic field-effect- 
transistor electronic systems has been hampered by the 
absence of simple and reliable models for their electronic 
quasiparticles and for the interactions of these quasipar- 
ticles with each other and with vibrations of the host lat- 
tice. Indeed the discovery of coherent band quasiparticlc 
properties and of the quantum Hall effect in these sys- 
tems, which are usually thought of as being complex and 
relatively disordered, has been one of the major surprises 
that has emerged from recent materials quality advances. 

The electronic properties of polyacene semiconductors 
are normally described by small polaron theory , which 
is quite successful at temperatures above «20 K. Small 
polaron theory starts from localized molecular orbitals 
and is therefore unable to describe the low-temperature 
hand quasiparticlc behavior seen in Battlogg et al. 's sam- 
ples, which have room-temperature to low temperature 
resistance ratios ~ 10'*. Schubnikov-de Haas transport 
studies in these samples establish quasiparticlc mean- 
free-paths in excess of 1000 lattice constants, demonstrat- 
ing that the low-temperature regime can be described 
only by starting at the opposite limit and building a the- 
ory based on delocalized band electron states. Never- 
theless, the interactions of quasiparticles with lattice vi- 
brations remain strong, and are presumably responsible 
for superconductivity. In this Letter, we propose a sim- 
ple tight-binding model for the quasiparticlc bands that 
emerge from 7r-molecular orbitals and are filled with car- 
riers by the field-effect, and for the interaction of these 



bands with the host crystal vibrational excitations. Out 
theory is parameter free once the band mass has been 
specified 

We find that over the entire wide range of studied 
densities, the quasiparticles of these systems lie within 
a single 2D band, bound to the crystal surface by the 
transistor's electric field. We argue that the dominant 
electron-phonon interactions in polyacene crystals differ 
qualitatively from the deformation and Frohlich interac- 
tions of inorganic semiconductors [Q and from the intra- 
molecule interactions that dominate in doped fuUcrenes 
being more similar instead to those of doped poly- 
acetelene |^ crystals. We describe these interactions us- 
ing a Su-Schrieffer-Heeger model detailed in the fol- 
lowing paragraphs, generalized to account for the larger 
number of important rigid molecule motions. The lattice 
phonons, a mixture in general of translational and libra- 
tional components, induce a change in the value of the 
hopping integral and hence couple to the quasiparticlc 
bands. We estimate the variation of hopping parameters 
with vibrational normal coordinates by assumming that 
they are proportional to the spatial overlap of Hiickel- 
type LUMO and HOMO molecular orbitals, with the 
proportionality constant fixed by matching the measured 
equilibrium hopping parameter [0. In this way, we have 
evaluated the carrier density dependence of a'^F{ijj), the 
index of electron-phonon interactions that figures promi- 
nently in the theory of phonon-mediated superconduc- 
tivity. We find peaks that match the spectrum extracted 
from recent tunneling experiments ||^. We also find that 
a^F((jj) is large, due to a density-of-states peak, when 
the lowest energy 2D subband is close to half-filling. We 
propose this property as the explanation of the relatively 
sharp onset of superconductivity at a 2D carrier density 
of w lO^^cm"^ seen in experiments [|l|. 

We focus in this study on anthracene which is mon- 
oclinic {P2i/a) with two basis molecules whose orienta- 
tions are related by a gliding plane symmetry |^. Since 
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Huckel model intra-molecular hopping parameters have 
no hnear dependence on low energy intra-molecular vi- 
bration normal coordinates, our model Hamiltonian in- 
cludes only inter- molecular hopping parameters B]. 



H = Tit + Tic-c 



Ti-v 



Ho 



where Ht is the LUMO or HOMO band Hamiltonian de- 
scribed by a nearest-neighbour hopping model, 7ic-o de- 
scribes Coulomb interaction between the electons, Tivib 
is the free-phonon Hamiltonian, and Tioxt describes the 
external electric field from the transistor's gate. 

We start by considering the field-induced energy bands 
obtained neglecting electron-phonon interactions. The 
key question we need to address is the number of 2D 
subbands that are occupied at a particular 2D density, 
n2D- As summarized in Fig. |], we find that all carri- 
ers reside in a single subband up to much higher density 
in these transistors than in their inorganic counterparts, 
principally because of the relatively small dielectric con- 
stants (e ~ 3.5) and the large effective in plane band 
mass {m± = l.Srrie) and in spite of the rather large c- 
direction effective band mass ~ 3me |p],p^ . Since the 
AI2O3 barrier is quite high (« 1.3 eV), we take it to be 
infinite. Our conclusion is based partially on envelope- 
function density-functional calculations for which the 
local-density-approximation (LDA) can be used to esti- 
mate many-body effects 0,11^ that favor a single sub- 
band. In our LDA calculations we find that the second 
subband is first populated at n2D ~ 5 — 7 x lO^^cm"^, 
a density just larger than the highest achieved in experi- 
mental systems. However, as shown in the inset of Fig. |^, 
the wavefunctions in these calculations are already local- 
ized over ^1 — 2 nm at densities well below this maxi- 
mum value. We have therefore repeated these electronic 
structure calculations using our tight-binding model, and 
treating interactions in a Hartree approximation. With 
this approach we find that 99% of the charge density is 
in the first molecular layer at n2D = 3 x lO^^cm"^, the 
density at which second subband occupation first occurs 
in this approximation. Since we believe that the exper- 
imental signatures of second subband occupation would 
be unambiguous and no anomalies in 2D density depen- 
dence have been reported, we conclude that a single 2D 
subband is occupied up to the highest densities and that 
the lowest subband is strongly localized in the top layer 
for n2D larger than « 10^'^ cm^^. 

To model the electron-phonon interaction, we expand 
each a — b plane hopping parameter to first order in the 
twelve coordinates that describe rigid rotations and dis- 
placements of neighboring polyacene molecules: 



m— 1 /i 



(1) 



where u^^m is the generalized displacement coordinate, 
II = 1,2 is the molecular basis index, m = 1,2,3 de- 



note displacements of the molecular center of mass along 
x,y and z directions, and m = 4, 5, 6 denote angular 
displacements around the 1,2,3 principal axes of each 
molecule. The electron-phonon interaction parameters 



dt/di 



are calculated based on the assump- 



tion of proportionality between hopping integrals and 
overlap integrals between HOMO or LUMO orbitals on 
adjacent molecules. The values obtained using standard 
Huckel approximation |Q HOMO tt orbitals (and hence 
appropriate to the hole systems on which we focus) are 
listed in Table I in units of io/a for m = 1, 2, 3 and for 
TO = 4, 5, 6. 
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SiOg/Si r\^^^^^^^=2. X lo'^cm ^ 
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FIG. 1. States per area with energy below the bottom of 
the second subband vs. 2D electron density for GaAs, Si, 
and polyacene field effect transistors from envelope-function 
LDA calculations. Second subband occupation occurs in 
the hatched area. The inset shows first and second sub- 
band wavefunctions obtained from envelope function LDA 
and tight-binding model Hartre-approximation calculations 
(bars) at n2D = 3 x lO^^cm"^. 

We now derive expressions for the interactions between 
phonons of the host crystal and the 2D band quasiparti- 
cles whose Hamiltonian is: 
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Qk6BZ2D j=l,2 



with 
e.(k) 
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cos( ^-^)-fcos(— 2 
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where Ea is obtained from the discrete layer subband cal- 
culation, the two electronic band indices originate from 
the two molecules per unit cell, and BZ2D indicates the 
two dimensional rectangular lattice Brillouin zone. The 
generalized displacements u^.m in Eq. ^ can be written in 
terms of quantized vibrational normal modes as follows: 



u,, 



.(R) = Ec 



X 
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Here N is the total number of unit cells, M^^m is the 
molecular mass of the /Ltth molecule for m = 1,2,3 and 
the molecular moment of inertia around the (m — 3)'th 
principal axis -f(^_3) for m=:4,5,6. The phonon frequen- 
cies and polarization vectors are obtained by solving the 
standard secular equation 



(k,q,g^,i/) = (-1) 2^ 



(4) 



These eigenmodes are in general a mixture of displace- 
ments and rotations [ p^ . We compute the dynami- 
cal matrix following the procedure outlined in Ref. p4[ . 
The phonon density-of-states and dispersion curves that 
emerge from these calculations are shown in Fig. ^j. The 
rigid-molecule vibration restriction used here is a conve- 
nient but inessential approximation that can be comple- 
mented if necessary by including coupling to important 
isolated molecule vibrations; it is however reasonably ac- 
curate for low-frequency vibrations of anthracene becom- 
ing less reliable for larger polyacenes. 
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FIG. 2. Calculated dispersion curves in anthracene (struc- 
ture shown in inset). Shown on the right is the phonon density 
of states FiLo) in units of meV~^ and normalized to 12, the 
number of vibrational modes considered. 

Combining these ingredients we finally rewrite the 
Hamiltonian (treating Coulomb interactions at the mean 
field level) as 



^ — ^2D + '^2Dc-vib + Wvib 

where 

'W2De-vib — 
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and 



xi^^iSy^M i.)e^q-^^.(e-»^"-[k+<il + (-l)(-^)e»*-'-k). (7) 

In Eq. ^/^^^(g^) = J2z, af*a'^e"i-'" is the form factor of 
the subband tight-binding wavefunction and the brackets 
indicate reduction to the 2D BZ. In our calculations we 
use the strict 2D limit {fa,a{qz) — 1) because of the 
strong quantum confinement found in LDA and Hartree 
calculations. 

We define the electron-phonon interaction spectral 
function in the standard way: 

a^F{oj)^j^ J2 5]ki;M(k,[k'-k],9.,^)P 

X ,5(e,(k) - ti)S{e,{k') ~ fi)6iLo - c^,([k' - k], g,) (8) 

where g{^) = J^k^i^iO^) ~ m) the electronic den- 
sity of states, and [k' — k] denotes the projection of 
the Q = (k' — k, g^) vector into the three dimen- 
sional BZ. The results of these calculations for densi- 
ties n2d = 4.0 X lO^'^cm"^, n2B = 1.4 x lO^'^cm"^, and 
?i2D = 6.7 X lO^'^ cm^-^ are shown in Fig. The peak 
locations in the low-frequency rigid molecule vibration 
regime agree with those observed in infrared absorption 
and tunneling data in pentacene crystals The mass 
enhancement factor A, obtained from the integration of 
2Fa'^{uj)/Lu, close to half filling is 0.25. Its sharp decay 
away from half filling (driven by the density of states) 
helps explain the sharp superconductivity onset observed 
in the experiments. We estimate the supercunducting 
critical temperature by the [dimension independent] BCS 
expression ksTc ~ Tiw/) exp[— 1/A]. Such expression, us- 
ing hcjo/kB ~ 150K from Ref. 1^, yields ^ 3K in agree- 
ment with the experiments We note that the inclu- 
sion of the second-nearest-neighbour hopping tends to 
push the density of states peak to higher energies (den- 
sities). Adding a finite life time of the quasiparticles, 
will broaden these peaks, hence creating a plateau in the 
n2D-dependence of Tc shifted to the right of half filling 
as observed in the experiments ||l[ . This last effect would 
lower also the Tc calculated at n = 4.00 x 10~^^cm~^, 
worsening experimental agreement. However, given that 
Ep does not satisfy Ep » TiLOmax, such estimates of 
from A are qualitative. 

In summary we have presented a theory of the low- 
temperature quasiparticle bands in polyacene field effect 
transistors, and of the coupling of these bands to the host 
molecular crystal phonons. We predict that the quasi- 
particles lie in a single 2D tight-binding band up to the 
highest densities that have been achieved experimentally 
at present and that the most important electron-phonon 
interactions arise from the influence of approximately 
rigid molecular translations and rotations on hopping be- 
tween HOMO and LUMO orbitals on adjacent molecules. 
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This picture implies a dependence of superconducting 
critical temperature on molecular lattice constant which 
contrasts with the case of doped fuUerene superconduc- 
tors. For the fuUerenes, the important electron-phonon 
interactions are intramolecular so that Tc depends on in- 
termolecular hopping only through the density-of-states. 
Decreasing the lattice constant, increases hopping, and 
hence decreases the density-of-states and Here 
decreasing the lattice constant will also strengthen the 
important electron-phonon interactions. According to 
our theory, the latter effect dominates and Tc will in- 
crease with decreasing lattice constant. We find that 
phonon-mediated electron-electron interactions in poly- 
acene molecular crystals are strong only when the 2D 
tight-binding band is close to half-filling and its density 
of states is relatively large. This occurs for 2D densities 
comparable to 4 x lO^^cm"^, where superconductivity 
turns on relatively abruptly. At lower carrier densities, 
A ^ 1 and electron-electron interactions are dominated 
by repulsive Coulomb interactions, consistent with the 
occurrence of the fractional quantum Hall effect. 
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FIG. 3. a^F{uj) 
n2d = 1.4x10^'' cm 
(dashed line). 



for n2d = 4.0 x 10^" c 
(dotted line), and n2d 



I ^ (solid line), 
6.7x10" cm"^ 
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TABLE I. Electron-phonon interaction parameters. Note 
that interaction parameters for different near neighbors are 
different. All interaction parameters are related by symmetry 
to the nine values list above. A complete list of hopping 
parameters for all near neighbors of the two molecules in the 
crystal's primitive cell is available from the authors. 
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